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The supercooled dynamics of a Lennard- Jones model liquid is numerically investigated studying 
relevant points of the potential energy surface, i.e. the minima of the square gradient of total 
potential energy V . The main findings are: (i) the number of negative curvatures n of these 
sampled points appears to extrapolate to zero at the mode coupling critical temperature Tc\ (li) 
the temperature behavior of n{T) has a close relationship with the temperature behavior of the 
diffusivity; (Hi) the potential energy landscape shows an high regularity in the distances among the 
relevant points and in their energy location. Finally we discuss a model of the landscape, previously 
introduced by Madan and Keyes [J. Chem. Phys. 98, 3342 (1993)], able to reproduce the previous 
findings. 
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INTRODUCTION 



The supercooled liquid regime is the interesting state 
of the matter which precedes the glass formation (once 
the crystal state is avoided). Many efforts have been 
devoted to better understand the relevant physical pro- 
cesses taking place in the supercooled state, the most im- 
portant being the enormous increase of both relaxation 
times and inverse diffusivity by many order of magnitude 
upon decreasing the temperature |l[] . In recent years the 
numerical investigation of simple model liquids has pro- 
vided a very useful microscopic description of the super- 
cooled regime and a detailed picture of the structural 
rearrangement of atoms during the dynamic evolution. 

One of the most powerful frameworks for the study of 
supercooled state has been the so called Potential Energy 
Surface (PES) description of the system [ 2|, ^ . This ap- 
proach focuses on the properties of the multi-dimensional 
surface of the total potential energy sampled during the 
time evolution of the representative point in the 3iV- 
dimensional configuration space (TV is the total number 
of particles). Different landscape features may control 
the motion of the particles @, | |, @, § |, 0, 
and the aim is to find those characteristics of the PES 
that have a direct relationship with (if possible the main 
causes of) the emergent behavior of the relevant physi- 
cal quantities. The challenge is to find the "good^ land- 
scape features that one supposes to be responsible for 
the interesting phenomena, or at least useful to give a 
clear interpretation of them. Among others, two differ- 
ent landscape approaches have been widely used in the 
last decades in this context: the Instantaneous Normal 
Modes (INM) approach based on the investigation 
of the PES very close to the instantaneous point in the 
configuration space during the molecular dynamics evo- 
lution of the system, and the Inherent Structures (IS) 
approach based on the analysis of the minima of 
PES visited by the system during its dynamic evolu- 



tion. In this paper we describe and revisit in a detailed 
way an additional approach, that has been very useful to 
give new insight in the analysis of the relevant processes 
taking place in the supercooled liquid regime. This ap- 
proach focuses on the minima of the square gradient of 
the total potential energy, "closest" to the instantaneous 
points of the molecular dynamics trajectory. This ap- 
proach allows one to obtain a microscopic interpretation 
of the relevant processes, the main result being the char- 
acterization of the dynamics above and below the Mode 
Coupling Temperature Tc iQ. Moreover, the analysis of 
these points allows one to obtain information about some 
relevant characteristics of the PES, of great importance 
to construct simplified models of the landscape. 

Lets us briefiy resume the two main PES approaches, 
the INM and th IS approach. The INM method is based 
on the investigation of the PES around the instantaneous 
configurations r (r represents the SA'^-dimesional vector 
of the representative point in the configuration space) 
during the molecular dynamics evolution of the system. 
The diffusive quantities are supposed to be related to 
the shape of the energy surface at r, that is to say to 
the eigenvalues and eigenvectors of the Hessian matrix 
(the second derivative of the potential energy). The 
main hypothesis of the INM approach is that the rele- 
vant diffusive directions have to be searched among the 
downward curvatures (eigenvectors with negative eigen- 
values). Many attempts have been devoted to extend 
this approach to different liquid systems and to develop a 
theory of the supercooled liquid state based on INM con- 
cepts. Moreover many efforts have been spent to recog- 
nize the true diffusive directions among all the downward 
ones, as it is known there are downward curvatures that 
do not correspond to diffusive directions, notably in the 
crystalline state. The INM negative curvature directions 
are classified as shoulder modes (related to anharmonic- 
ities of the PES) and double wells (with a double well 
shaped one-dimensional profile) . Diffusive directions are 
finally identified as those leading to different minima pq] 
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(this analysis involves the search of the minima of the 
PES). The temperature dependence of the diffusion coef- 
ficient has been shown to be related to the temperature 
dependence of the fraction of diffusive directions jl^ . On 
the basis of this analogy, it has been conjectured a struc- 
tural interpretation of the Mode Coupling Temperature 
Tc as the temperature at which the number of diffusive 



directions goes to zero |17 



The second IS landscape approach is based on the anal- 
ysis of the inherent structures, the minima of the PES. 
To each instantaneous configuration r one associates a 
inherent structure r^^ (using, for example, a steepest de- 
scent path starting from r), r — > r^^, partitioning the 
whole PES in the basins of attraction of minima. This 
kind of analysis has been very useful not only in the study 
of the diffusive directions in the supercooled regime (see 
the previous discussion about the INM), but also in the 
study of thermodynamic quantities (for example to evalu- 
ate the configurational entropy ^1) and in the study 
of the out-of-equilibrium dynamics 1|2^, . 

In the present work we adopt a new approach in the 
analysis of the PES, in some way intermediate between 
the previous two, in that it maintains information about 
the possible diffusive degrees of freedom (as the INM ap- 
proach) and uses a mapping in the PES that associates 
instantaneous configurations r to new points in the PES 
(as the IS approach). In this approach |2^, ^ the min- 
ima of the square gradient of the potential energy 
W — jVFp (saddles and some inflection points of the 
PES), reached starting from equilibrated configurations 
r, are calculated. This procedure partitions the PES in 
basins of attraction of this new points: r ^ r^. The 
choice of W, suggested long time ago by Weber and Still- 
inger , was motivated by the fact that saddles points 
of V are absolute minima of W . As pointed out recently 
by Doye and Wales |^ and discussed in more details in 
the following, the numerical minimization of W locates 
mostly local minima of W ^ i.e. points which are not real 
saddles of V . For this reason in the following we refer 
to the local minima of W as Quasi Saddle Points (QSP) 
pGf . In Rcf.s and the topological and metrical 
properties of such points have been used for a description 
of the dynamics of supercooled liquids. The investigated 
systems were monatomic Lennard- Jones and/or Binary 
Mixture Lennard-Jones. After equilibration of the sys- 
tem at a given temperature T, corresponding to a given 
instantaneous potential energy e, the minima of W were 
searched and associated to the stationary points of V 
(saddles) . Then the energy and the number of nega- 
tive curvatures at these saddle points were measured. 
The dependence of and on the temperature ||2^ 
or on e | |2^ show interesting properties: (i) the quan- 
tity extrapolates to zero at Tc, demonstrating the 
validity of the conjecture that Tc marks the transition 
between a dynamics among minima at low T and a dy- 
namics among saddles at high T; {it) the aspects of the 
energy landscape seen by a given minimum is highly reg- 
ular (as demonstrated by the linear dependence of on 



— or on [|3|)- "^^^ previous observation 

leads to the conclusion that the knowledge of the saddle 
points properties can be used to predict the supercooled 
dynamics and points out the relevance of these peculiar 
PES points. 

Following these papers, the concept of saddles - as ob- 
tained by a minimization of |Vl^p - has been further 
used in the investigation of the supercooled liquid dy- 
namics ^ and aging dynamics |^ . 

As we anticipated above, there is a drawback in the 
approach used in Ref.s and ||2^. In fact, while it is 
true that in a saddle point (stationary point) of the 3A^- 
dimensional surface V{v), the function W has a minimum 
(actually, W — there and W has a global minimum), 
the reverse is not guaranteed to be true. There are lo- 
cal minima of VF, with W 0, that are NOT stationary 
points of V. In these points, along one or more directions 
r(^), parametrized by ^, the function V{^) has an inflec- 
tion point. In these inflection points (QSP), the Hessian 
of V has a number of vanishing eigenvalues (excluding 
the three translation directions) equal to the number of 
directions r(^) where V has an inflection point. This 
property has been used in in order to identify the 
QSP: a threshold was arbitrary chosen, and all the mini- 
mizations of W that lead to points where the Hessian of 
V has at least one eigenvalue smaller (in modulus) than 
the threshold were discarded from the analysis. In 
we claimed that only a small fraction of minimizations 
ended in a QSP points, and we estimated in 2% the error 
introduced in by the choice of the threshold. 

Doye and Wales ||2^ have shown - by the analysis of the 
same system investigated in 22 and , the BML J with 
periodic boundary conditions - that almost all (around 
95% in their work) the minimizations of W ended in a 
QSP point, and not in a true saddle as claimed in 
and |23). Motivated by such contradictory results, we 
reanalyzed the minima of W found in . In this work 
we show the result of such an analysis. We conclude that 
the results of Doye and Wales are correct and the error in 
the analysis in originates from the "imperfect" min- 
imization of W: in a point close to (but not exactly at) 
QSP, the Hessian of V no longer has "small" eigenvalues, 
and the criterion used to identify false saddles fails. 

One puzzling question, however, still remains open: 
why the quantity Ug - i.e. the number of negative eigen- 
values at the minima of - is related to the critical tem- 
perature of the system ? It remains true that (T) ex- 
trapolates to zero at Tc- In the present paper we address 
this question, and we will further stress the importance 
of the quantity ; we will show that also the transport 
properties (namely the diffusion coefficient) can be de- 
termined by the knowledge of . 

The paper is organized as follow: In section II we de- 
scribe the new approach, defining the saddles of the PES 
and studying very carefully their operative definition and 
the problems related to the presence of false saddles. 
Then in section HI we expose the relevant results ob- 
tained with the new PES approach, giving a topological 
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interpretation of the mode coupling dynamic transition. 
In section IV the relationship between difFusivity and sad- 
dle properties is evidenced, giving a simple interpretation 
of the diffusion processes in terms of saddle order. In sec- 
tion V we describe some landscape features as emerged 
from saddle analysis. In section VI a simplified model 
of the PES is described, together with its relevant prop- 
erties compared to the features of model liquids investi- 
gated above. At the end (section VII) a brief summary 
of the main results and conclusions are reported. 



II. QUASI SADDLE POINTS 

The analyzed system is a simple model liquid, a 
monatomic Modified Lennard-Jones 6-12 (MLJ). The 
model is able to support strong supercooling without the 
occurrence of crystallization, due to the presence of a 
small perturbation term in the Hamiltonian that inhibits 
ordering (this term is a function of the static structure 
factor, see Ref. for more details). The system is com- 
posed of N — 256 particles enclosed in a cubic box with 
periodic boundary conditions. Truncated (i?c — 2.6) and 
shifted LJ potential are used. Equilibrium configuration 
are prepared by standard microcanonical molecular dy- 
namics simulations at constant density p = 1 (standard 
LJ units are used hereafter) and at temperatures rang- 
ing from the normal liquid phase (T ^ 1.6) down to the 
mode coupling temperature {Tc = 0.475, estimated from 
diffusivity) . 

The aim of our approach is to use a method of inves- 
tigation of the PES that allows us to reject all the non- 
relevant degrees of freedom for the description of the slow 
supercooled dynamics. The idea is that there are partic- 
ular points of the PES associated in some way to the 
instantaneous configurations, and that the properties of 
these points contain all the important information for the 
determination of the slow dynamics. Possible and natural 
candidates are the saddle points of the PES. In principle 
one can think to partition the whole configuration space 
in basins of attraction of saddles (the definition of the 
partition is obviously not unique, a possible choice could 
be that obtained using the Voronoi polyhedra in the 3iV- 
dimensional space) . However a useful partition in basins 
of saddles (that is local, i.e. allows us to associate a sad- 
dle point to each instantaneous configuration using only 
local information of the PES) is still not available. A pos- 
sible partition is that obtained using the pseudo-potential 
Q W = |VT^p, associating to each instantaneous con- 
figuration during the molecular dynamics a minimum of 
W: 

instantaneous configuration > minimum of W . 

However, as already pointed out, the relation between 
all the minima of W and saddles of V is not one-to-one: 
absolute minima of W (with W = 0) are all saddles of V 

absolute minima of W < > saddles of V , 



but local minima (with > 0) are inflection points of 
V (more precisely they are inflection points with the first 
derivative of V of the same sign of the first derivative of 
the curvature in the inflection direction): 

local minima of W > inflection points of V . 

In Fig.s |l| and ^ examples of two possible one-dimensional 
profiles of V and corresponding W are shown as a func- 
tion of a generic one-dimensional coordinate. In Fig. ^ 
the two absolute minima of W correspond to two saddle 
points of V (in one dimension a minimum and a max- 
imum). The arrows define the basins of attraction of 
minima of W, and the boundaries of these basins are in- 
dicated with open symbols. In Fig. ^jthe local minimum 
of W corresponds to an inflection point of V. 
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FIG. 1: Example of profile of V and W along a given direction. 
The minimum and maximum of V correspond to absolute 
minima of W . The arrows indicate the basins of attraction of 
the minima of W . 

We note that at inflection points the eigenvalues of the 
Hessian vanish, so in principle one can recognize them 
simply analyzing the eigenvalues of the Hessian of V . 
In a previous work ]2^ we used this eigenvalues crite- 
rion (choosing a threshold to discriminate inflection di- 
rections) and we found that local minima of W were 
found very rarely. A closer inspection, however, has re- 
vealed that, due to "imperfect" minimizations of W ^ the 
eigenvalues associated to a point "close" (due to the not 
perfect minimization) to an inflection point were of the 
same order of magnitude of the eigenvalues associated 
to a saddle. Refining the analysis of the minima of PF, 
we find that the number of points that are true saddles 
is very small (about 1% of the total), so almost all the 
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FIG. 2: Example of profile of V and W along a given direction. 
The inflection point of V corresponds to a local minimum of 
W with W > 0. The arrows indicate the basins of attraction 
of the minima of W . 



minima of W found are inflection points. However there 
is another interesting characteristic that has to be con- 
sidered, the number of inflection directions (number of 
zero eigenvalues of Hessian, excluding the tranlations) at 
the minima of W: if this number is small one can think of 
the minima of as a "quasi" saddle point QSP, as they 
are true saddle points in the subspace orthogonal to the 
small number of inflection directions. One can conjecture 
that the properties of a QSP are very similar to that of 
a true saddle point associate to an instantaneous config- 
uration; however, a relationship between them has still 
to be proved. In the simulated MLJ system the number 
of inflection directions uq with respect to the number of 
negative curvatures are reported in Table || for the dif- 
ferent investigated temperatures. The values of tiq (from 
I to 4) indicate that they correspond to few directions 
in the configuration space; a small value that allows us 
to speak about QSP as good candidates to approximate 
the properties of the true saddle points. We note that 
the fraction / = no/n^ is higher for the low temperature 
data (where the quantity is small), evidencing how 
these points are more influenced by inflection directions. 
However, also in this case it is possible to estimate the 
error in the calculation of the order (see Fig. || in the 
next section), evidencing the robustness of the obtained 
results. 

A natural question now arises: what are the relevant 
features of the found minima of W that are useful to de- 
scribe the long time dynamics of the supercooled liquids 



TABLE L Average number of negative curvatures (n^) and 
inflection directions (no), and their ratio (f) at minima of W 
for different temperatures. 



T 


n" 


no' 




.49 


1.4 


1.1 


0.78 


.53 


3.0 


2.4 


0.80 


.57 


5.6 


2.5 


0.45 


.61 


7.0 


3.1 


0.44 


.64 


7.6 


2.8 


0.37 


.71 


10.3 


2.9 


0.28 


.81 


15.9 


3.6 


0.23 


.97 


20.3 


3.4 


0.17 


.13 


24.5 


4.3 


0.18 


.28 


30.1 


3.9 


0.13 


.44 


32.9 


3.9 


0.12 


.57 


36.2 


4.4 


0.12 



"Absolute number of negative curvatures 

''Absolute number of zero curvatures (excluding translations) 
'^Fraction of inflection directions no/n^ 

analyzed ? In the spirit of our early simple conjecture 
that the relevant quantity is related to the number of 
negative curvatures at saddle points, we speculate that 
the relevant information is in the negative curvatures and 
not in the few inflection directions. So, no matter if they 
are true or quasi saddle points, we are going to check if 
their properties are able to describe and to give a mi- 
croscopic interpretation of the dynamic processes in the 
supercooled regime. 



III. DYNAMIC TRANSITION 

We firstly analyze the temperature behavior of some 
quantities related to quasi saddle points, and compare 
them to the corresponding behavior of the same quan- 
tities calculated at instantaneous configurations and at 
inherent structures. 

For each temperature (ranging from T ~ 1.6 to T ~ 
Tc), we analyze 20 independent equilibrium configura- 
tions. For each configuration we calculate the associ- 
ated IS, implementing a steepest descent algorithm which 
moves in the direction of —WV = F, and the associated 
minimum of W, moving in the direction of —\/W = 7i-F, 
where Ti. is the Hessian matrix. So, at each temperature, 
we have three different points of the configuration space: 
the instantaneous configuration r, the inherent structure 
r^-g, and the quasi saddle configuration r^. 

In Fig. ^the temperature dependence of the energies of 
the different points are shown: e, e^g and e^. The value 
of the mode coupling temperature Tc = 0.475 (obtained 
by diffusivity data, see section IV) is also indicated. The 
energy of the inherent structures shows a constant behav- 
ior down to T ~ 0.8, below which it starts to decrease 
abruptly (note however the small variation range with 
respect to that of the e^,). Below this temperature the 
system starts to visit basins of minima of lower and lower 
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FIG. 3: Instantaneous e, saddle , and inherent structure 
Cjg energies as a function of temperature. 



energy. The quantity shows a different behavior, de- 
creasing gradually with temperature and approaching e^^ 
close to Tc- The fact that lies well below e indicates 
that the process of minimization of the pseudo-potential 
W consists in a downward path in the PES along the 
majority of the degrees of freedom. The coincidence of 
e^g and at Tc suggests the hypothesis that the sam- 
pled saddles start to be mainly minima at this temper- 
ature. The answer to this question is in the behavior of 
the number of negative curvatures as a function of tem- 
perature. In Fig. ^ the fraction of negative curvatures 
(number of negative curvatures over the total degrees of 
freedom 3iV) is shown for the analyzed configurations r 
(n) and (n^) (obviously this number is zero for the 
inherent structures r^-g). The first thing to observe is 
that the quantity lies well below the instantaneous 



order n, indicating that the process of minimization of 
W leads to a point well below in the PES, with a number 
of negative curvatures less than that at the equilibrium 
starting point. We have then found a direct method to 
reject many negative eigenvalues (corresponding to not 
diffusive directions) of the Hessian at the instantaneous 
configurations. If the remaining negative curvature di- 
rections are those useful to describe the slow dynamics 
of the system, will be soon clear. A first good indication 
of the relevance of the quantity for the description of 
the dynamics is the fact that its well defined temperature 
behavior extrapolates to zero very close to the Tc value. 
In Fig. ^ the temperature behavior of (T) is shown 
in a larger scale. It is shown the calculated in a di- 
rect way (as the fraction number of negative curvatures 
at the minima of W) and an upper and lower estimation 
of the "true" saddle order, considering all the inflection 
directions as if they contributed to the saddle order (up- 
per points) and as if they didn't contribute to it (lower 
points). Fitting the three different data with a power 
law (T — To)''', we obtain very similar value of Tq ~ Tc, 
indicating the robustness of the analysis. 



0.06 





FIG. 4: Instantaneous order n, and saddle order as a 
function of temperature T. 



FIG. 5: Saddle order function of temperature. The 

open squares are the upper and lower estimations considering 
the "false" directions (see text for details). 

The sampled minima of W are then QSP with a num- 
ber of negative curvatures that decreases lowering the 
temperature, and approaches zero at the temperature of 
the dynamic mode coupling transition. At Tc and below 
it, T < Tc, the sampled saddles are in average minima of 
V, i.e. saddles with order zero. The emergent scenario 
for the dynamics is then the following: 

• for T > Tc the system lies with high probability 
close to borders of basins of attraction of inherent 
structures, 

• for T < Tc the system spends most of the time 
trapped in the basins of inherent structures. 

These findings confirm the conjecture that Tc marks a 
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dynamics crossover, evidencing a possible structural PES 
interpretation of this change in the dynamics. 

We conclude this section showing the distributions of 
the quasi saddle order for different temperatures. In Fig. 
^ the distribution P„ of saddle order n (now we indicate 
with n the absolute saddle order, i.e. the total num- 
ber of negative curvatures, not the fraction) is shown for 
four different temperatures. A simple conjecture allow us 
to obtain a prediction for these distributions. Assuming 
the independence of the relevant slow degrees of freedom 
and a simple one dimensional profile of the independent 
coordinate (with only a minimum and a maximum) we 
have that the distribution of the saddle order n of N 
degrees of freedom is simply the binomial distribution 
Pn = (^)p"(l — where p is the probability that 

the single degree of freedom is close to (inside the basin 
of attraction of) the maximum. In the large N and small 
p limit we obtain a Poisson distribution 



< n >" 

r(n + i) 



-<ri> 



(1) 



written in term of the P function (P(n + 1) = n\) and 
of the mean value of the saddle order < n >— p N (the 
temperature dependence is in the p parameter - small p 
means low temperature). Lines in Fig. ^ are the distri- 
butions obtained from Eq. |^, with the < n > value fixed 
by the data. They well describe the measured distribu- 
tions, suggesting the correctness of the simple viewpoint 
of independent degrees of freedom (a simple model with 
this ingredient will be study in detail in the last section). 



system are able to capture the dynamic crossover at T^. 
In order to corroborate this finding we now go to analyze 
the diffusion properties in the supercooled regime and to 
evidence a possible relationship with the saddle order . 

For different temperature above Tc we have calculated 
the mean square displacement: 



> 



(2) 



and, from its long time behavior, the diffusion coefficient 
D: 



(3) 



In Fig. the diffusion coefRcient D is shown as a function 
of temperature. The mode coupling theory predicts a 
power law divergence of the inverse diffusivity close to 
Tc, so the behavior of the diffusion coefRcient is a method 
to estimate Tc- Fitting the data with power law (line in 
Fig. I) 



D{T) cx (T - T,) ' 



(4) 



we obtain the value Tc — 0.475. 

The fact that also the saddle order extrapolates to 
zero at Tc, suggests the possibility of a relationship be- 
tween and the diffusivity D. A direct prediction of a 
relation between the two quantities arises from a simple 
interpretation of the dynamic processes. Supposing that 




n 



FIG. 6: Distribution of saddle order n (the absolute number 
of negative curvatures) for four different temperatures. The 
line are the Poisson distributions (see the text). 



IV. DIFFUSION 

As demonstrated in the previous section, some features 
of the saddles probed during the dynamic evolution of the 
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FIG. 7: Diffusion coefficient D as a function of temperature. 
The line is a power law fit with parameters reported in the 
text. 

the slow degrees of freedom responsible of the diffusion 
processes are in some way related to the number of neg- 
ative curvatures of the probed saddle points at a given 
temperature, one can try to predict a relationship be- 
tween the saddle order and the diffusion coefRcient. 
Assuming that the number of diffusive directions at a 
given temperature are proportional to the saddle order 
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at the same temperature, we can think of the diffusive 
process as a random walk process in a space of dimen- 
sion equal (or proportional) the the saddle order . In 
this case one expects that the diffusion coefHcient is pro- 
portional to the temperature times the saddle order: 



D{T)^Tn,{T) 



(5) 



We note that this hypothesis doesn't mean that in gen- 
eral all the negative curvatures at a generic saddle point 
are related to a slow diffusive process, as there are sad- 
dles related to small local rearrangements of atoms. It 
means only that the process of minimization of VF, start- 
ing from an equilibrated configuration, leads in average 
to a point (saddles or quasi-saddle doesn't matter) with 
negative curvature directions in some way related to the 
true diffusive directions. 



V. LANDSCAPE FEATURES 

The previous analysis of the PES allows us to infer 
some topological characteristics of the landscape itself. 
The first feature of the PES we analyze is the energy or- 
ganization of the saddles, i.e. a possible relation between 
energy and order of saddle points. Plotting (see Fig. ^a) 
the saddle energy eg as a function of the saddle order 
(using all the saddles found at different temperatures), 
i.e. 



as done in |23 



;("s) « ^e' rig 



(6) 



the data can be reproduced with Ae' = 13.5. However, 
data at low energy (and order) deviate from the linear 
fit. Observing that these points are those sampled in the 
temperature range where the inherent structures start to 
decrease in their energy values, one is attempted to re- 
late the deviation from the linear fit of the energy-order 
relation to the inherent structure behavior. In Fig. ^ 




T 



FIG. 8: Diffusivity over temperature D/T and saddle order 
rig (multiplied by a constant factor a) as a function of tem- 
perature. The line is the power law fit of diffusivity. In the 
inset the two quantities one against the other in a double log 
scale. The line has slope equal to 1. 



In Fig. 1^ the diffusivity over the temperature calcu- 
lated through molecular dynamics simulations and the 
saddle order are shown as a function of temperature: 
it seems that the simple random walk model on a reduced 
space works well. Eliminating the temperature and plot- 
ting the two quantities one against the other in a double 
log scale (inset of Fig. ^) we see that the Eq. || is well 
suited, except the last low temperature point, which is 
however the most influenced by a possible not good equi- 
libration procedure. This kind of analysis deserves surely 
further investigations, for example analyzing in a more 
direct way the paths during the molecular dynamics evo- 
lution and their relation with the negative curvature di- 
rections at saddle points. However one result emerges in 
a strong way: the quantity seems to be very useful to 
reproduce dynamic processes and to give a possible PES 
interpretation of them. 
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FIG. 9: Saddle energy (a) and saddle energy elevation 
from underlying minima — e^g (b) as a function of saddle 
order . 

the quantity plotted against the saddle order is the 
elevation energy — e^^g of saddles with respect to the 
corresponding local minima (minima visited at the same 
temperatures), obtaining now a remarkable linear rela- 
tionship in the full temperature range: 



Ae 



(7) 



with Ae = 13.3. This suggests that the energy landscape 
above a given minimum is organized in families of equally 
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spaced energy saddle points, with only a single energy 
barrier parameter (in average) Ae, that represents the 
energy gap between a saddle of order and a saddle 
of order + 1. A possible explanation of this result 
is the following: saddles of a given order are obtained 
from a combination of independent saddles of order 1, 
that lie a fixed amount over minima. One can think of 
the excitations as as a "gas" of non-interacting degrees 
of freedom (in the next section we develop this viewpoint 
in a deeper way, introducing a simple landscape model). 
Similar results are found in the out-of-equilibrium regime 
pof , suggesting that the PES properties obtained are not 
influenced by the kind of dynamics the system use to 
explore its PES. 

Another interesting topological information about the 
PES is obtained calculating the distance relations be- 
tween adjacent inherent structures. For each saddle point 
we have perturbed the system along a randomly chosen 
negative curvature direction and then started a mini- 
mization procedure of the potential V in order to find the 
underlying minimum. The same minimization has been 
performed perturbing the system along the previous neg- 
ative curvature direction but in the opposite versus. In 
this way we obtain couples of inherent structures asso- 
ciated to each negative direction (we call them adjacent 
inherent structures). We then evaluated the Euclidean 
distance between adjacent IS, defined as 



1 



dt. = — 



N ^ 



< \n{ISi)-n{IS2)\^ > 



(8) 



where ISi and IS2 are two adjacent minima and the 
average is over different couples of adjacent IS. In Fig. 
[l0| the quantity d^^ is shown as a function of the order 
rig of the starting saddle. The observed linear relation 
between the two quantities 



(with dg = 9.6) suggests a simple topological interpreta- 
tion: the descent path from a saddle of order towards 
the underlying minima can be view as a sequence of in- 
dependent random steps, each of them decreasing the 
saddle order by 1 and increasing the mean square dis- 
tance between the underlying local minima by a constant 
amount dp. 

The above results seem to indicate that in some aspects 
the PES exhibits a very simple and organized structure. 
This suggests the possible use of simple solvable mathe- 
matical models in order to capture some of the relevant 
features of the PES explored in the supercooled regime. 



VI. TRIGONOMETRIC MODEL 

A very simplified model that is able to capture some of 
the regular characteristics of the PES of liquids evidenced 
above, is the so called Trigonometric Model (TM), intro- 
duced by Madan and Keyes js^ . In this section we calcu- 
late for the TM the behavior of the quantities previously 
analyzed for the Lennard- Jones system. 

The TM is a model for N independent degrees of free- 
dom with Hamiltonian: 



ff,,, =A^[l-cos(^0] 



(10) 



dUns) 



dn 



(9) 



where {ipi} are angular variables: (pi G [0, 27r). The PES 
of TM reproduce the regularity of the average saddle 
properties of the Lennard- Jones PES: saddles of order 
are 2A over (in energy) and tt near (in distance) saddles 
of order — 1. The thermodynamics is easy computable, 
as the partition function is factorized: 

Z,„(/3) = Zo~(/3) , (11) 

where (3 — (we use the unit Kb = 1) and 
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FIG. 10: The mean square distance between adjacent minima 
vs. saddle order n^. 



dip e-'^^Ii-^^^^)! = 2n e'^^ /o(/3A) , 

(12) 



where Iq is the Bessel function of order zero. The energy 
density e^^j^ = —N^^dp log(Zyj,j) is easy written in term 
of Iq and Ii (the Bessel function of order 1, Ii{x) = 



1 



lom) 



(13) 



In order to test the reliability of TM in the description 
of liquid behavior, we now calculate the diffusion coeffi- 
cient and the saddle order as a function of temperature. 
The diffusion coefficient D^^^^ is exactly computable cal- 
culating the mobility in linear response | |3^ , with the 
introduction of a Langevin dynamics: 



(14) 
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FIG. 11: The reduced diffusivity D^j^^/Dq as a function of re- 
duced temperature T/A. The symbols are from iso-potential 
simulation (see the text). 



where F is an external force, 7 is the friction, 77 is a 
random variable with mean zero and (5-correlate in time: 
< rj{t)ri{t') >= {2jjf3)5{t — t') (we choose m = 1 mass 
unit). Following p3] we obtain: 



1 



(15) 



Using the expansion of /g for large argument ( /o(a;) 
exp(a;) / ^/x, for large x) we obtain the Arrhenius behavior 
for the diffusivity at low temperature: 

D^^, cx e"'^'^^ for low T (16) 

(we note the value 2A is the energy barrier value for the 
single degree of freedom). The dimensionless quantity 



Do 



-'0 



(A/3) 



(17) 



where Dq = (7/?)"^ is the free Brownian diffusion co- 
efficient, is shown in Fig. |ll| as a function of temper- 
ature (line). An interesting observation arises from the 
comparison between this theoretical canonical result and 
a numerical iso-potential (dynamics at constant poten- 
tial energy) calculation. The points in Fig. |ll| are ob- 
tained from a simulation of N = 1000 variables con- 
strained to move at constant potential energy (the dy- 
namics is a random walk dynamics in which the forces 
— Vi-ff are projected to the constant potential energy sur- 
face) . The coincidence between analytical and numerical 
data evidences how canonical and iso-potential calcula- 
tions lead to the same equilibrium dynamic properties. 
In Fig. |l^ the diffusion coefficient D.j,m shown as a 
function of temperature in a semi-logarithmic scale. We 
note that, though D^^ is a smooth function, one can 
define a change in the behavior extrapolating a critical 
temperature Tc, using a power law fit in the low temper- 
ature range (from 0.5 to 1.5): 

D,,,{T)^{T-T,y , (18) 



10 



S 10" 



10 ■ 



A (T - T )^ 




0.0 



FIG. 12: The diffusivity -D^j,, as a function of reduced tem- 
perature T/A in a logarithmic plot (full line). The dashed 
hne is the power law fit of the high temperature range 
(Tc/A = 0.38, 7 = 1.2). In the inset the low temperature 
behavior of diffusivity -D^j^ is shown as a function of inverse 
temperature in an Arrhenius plot (dashed line indicates slope 
2). 



obtaining the values: Tc/A = 0.38 and 7 — 1.2. Below 
Tc the diffusion coefficient can be approximated by Eq. 
[1^ , a low T Arrhenius behavior with energy barrier value 
2 A (see the inset in Fig. 

The saddle order n^^j^ is also computable in an exact 
way. Defining it as the negative curvatures at minima of 
W = |Vi/p (we note that, in the TM, all the minima of 
W are true saddles of H and, moreover, the instantaneous 
order coincides with the saddle order) we can calculate 
the equilibrium mean value of rir^j^ as the probability 
that the variable ip is in the range [7r/2, 37r/2], where the 
H profile as a negative curvature: 



1 / dV? e"'''"'~'°''^^^' 



Tr/2 



f3iT/2 



The temperature behavior of rir^^^ is shown in Fig. [l^ . 
Also in this case it is possible to fit the data with a power 
law fit close to Tc (from 0.5 to 1.5), obtaining a value 
of Tc/A = 0.41, close to that obtained from diffusivity 
(Tc/A = 0.38). 

Having calculated the temperature behavior of diffu- 
sivity D^^j{(3) and saddle order n^^j{l3), we can now try 
to study the relationship between them. Due to the fact 
that a direct analytical expression of D^^^^ as a function of 
n^j,^ is not available, we study it numerically eliminating 
the temperature parameter. The relationship is shown in 
the inset of Fig. |3[ where the quantity D^^j/Dq is plot- 
ted against n^^j- in a logarithmic scale (this plot corre- 
sponds to that of the Lennard- Jones system - inset of Fig. 
- as the quantity Dr^^^j / Dq is proportional to f]D^j^,/). 



10 




FIG. 13: The saddle order a function of reduced 

temperature T/A. The dashed hne is the power law fit in 
the temperature range 0.5 - 1.5 (T^/A = 0.41, 7 = 0.58). 
In the inset reduced diffusivity D^^^/Do against the saddle 
order n^j^ in a double logarithmic scale. The dashed line has 
slope equal to 1. 



The investigated T range is close to the critical tempera- 
ture Tc (as in the liquid system) and also in this case we 
find an approximatively linear relationship (the dashed 
line in the inset of Fig. ^ corresponds to a power law 
with exponent 1): 



D„., (X T 



(20) 



In conclusion, the very simple TM is able to reproduce 
some of the properties emerged in the analysis of the 
Lennard-Jones liquid system, indicating that the use of 
simplified analytical model of the PES is a promising field 
of research. 



VII. CONCLUSIONS 

The analysis of the PES has been very useful in the 
investigation of supercooled dynamics of simple model 
liquids. Besides the two main approaches in the investi- 
gation of the PES, the instantaneous normal mode and 
inherent structures approaches, we have discussed an ad- 
ditional sampling of the PES explored during the dy- 
namic evolution of the system in configuration space: it 
consists in the analysis of the minima of the square gradi- 
ent of the potential energy, W = |Vl^p. The character- 
istics of these sampled points seem to have a very close 
relationship with the transport properties of the system. 
The main result we obtained is the characterization of the 
dynamic transition temperature Tc as the temperature at 
which the behavior of quantities related to the minima 
of W changes on character: more specifically the number 
of negative eigenvalues of the Hessian at these points ex- 
trapolates to zero at Tc. The minima of W sampled at 
different temperatures appear to contain the relevant in- 
formation about the diffusive directions and the process 



of minimization of W eliminates the non-diffusive neg- 
ative directions usually present at instantaneous config- 
urations. Recently, an interesting technique to evaluate 
the non-diffusive directions as been adopted in p3|] . In 
some previous works we called these points saddles, as 
the absolute minima of the function W are true saddle 



points of V. However, following Doye and Wales |25 



closer inspection has revealed that very often the sam- 
pled minima during dynamic evolution are local minima 
of W, that correspond to inflection points of V (the en- 
ergy profile along some normal mode directions at these 
points is an inflection direction). The fact that the num- 
ber of inflection directions at a given local minimum of 
W is small, allows us to think that the main properties 
of these points are determined by the saddle directions 
(that give the value of the order, i.e. the number of neg- 
ative curvatures). A further evidence of the relevance 
of the minima of W is obtained analyzing the diffusiv- 
ity. The temperature behavior of the number of negative 
curvatures is related to that of the diffusion coefficient, 
indicating a close relationship between them. It could be 
that the main information about the dynamic processes 
is in the saddle points and the relevant properties of the 
local minima of W are good approximations of the prop- 
erties of these true saddle points. But there is also the 
possibility that the true saddles are not so relevant in the 
description of the dynamics, at least not more relevant 
than points with inflection directions (local minima of 
W). 

All these findings lead to the following viewpoint about 
the supercooled dynamics: the relevant diffusion pro- 
cesses happen in a subspace of the total configuration 
space with dimension proportional to the number of sad- 
dle directions at the sampled minima of W, and, due to 
the presence of these "open" directions, the relevant dif- 
fusive processes are not activated in energy. Below Tc 
things are different, as the number of negative directions 
at sampled minima of W is vanishingly small, indicating 
that the system is with high probability in a basin of an 
inherent structure (minimum of V). In order to change 
basin and diffuse, now the system has to find the "good" 
directions leading to other minima (an entropic process, 
considering the whole SN dimensional space of the de- 
grees of freedom, that could correspond to an energy ac- 
tivated process in a suitable n dimensional subspace - 
with n < 3iV). 

Another important result obtained refers to the struc- 
tural organization of the PES, i.e. the relationship among 
the relevant points of it. A very simple structure of the 
PES emerges: saddles of order n -t- 1 lie above saddles of 
order n by a constant quantity Ae and adjacent inher- 
ent structures are equally spaced in configuration space. 
The above observed regularities suggest the possible use 
of simplified models of the PES in order to capture rel- 
evant characteristics of the supercooled dynamics of liq- 
uids. The model we analyzed here is the Trigonometric 
Model, introduced by Madan and Keyes a model of 
independent sinusoidal degrees of freedom. The behav- 
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ior of the calculated difFusivity and saddle order seems to 
indicate that the model, despite the absence of cooper- 
ativity, is able to reproduce important aspects of liquid 
dynamics close to Tc. The study of more complex and re- 
alistic PES models, with interactions among the degrees 
of freedom, is a interesting development of this kind of 
analysis. 
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